Ion-trap quantum computing in the presence of cooling 



Almut Beige 

Max- Planck- Institut fiir Quantenoptik, Hans-Kopfermann-Str. 1, 85748 Garching, Germany and 
Blackett Laboratory, Imperial College London, Prince Consort Road, London, SW7 2BW, UK* 

(Dated: February 1, 2008) 

This paper discusses ways to implement two-qubit gate operations for quantum computing with 
cold trapped ions within one step. The proposed scheme is widely robust against parameter fluctu- 
ations and its simplicity might help to increase the number of qubits in present experiments. Basic 
idea is to use the quantum Zeno effect originating from continuous measurements on a common 
vibrational mode to realise gate operations with very high fidelities. The gate success rate can, in 
principle, be arbitrary high but operation times comparable to other schemes can only be obtained 
by accepting success rates below 80%. 
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I. INTRODUCTION 

Very recently, considerable progress has been made 
in the experimental reahsation of quantum computing 
schemes with cold trapped ions. In Innsbruck, the Cirac- 
ZoUer controUed-NOT quantum gate Q has been imple- 
mented with the help of six concatenated laser pulses 
individually addressing each of the two ions 0. At the 
same time, the group in Boulder demonstrated a robust 
high-fidelity geometric two-qubit phase gate in the labo- 
ratory. This was achieved with a sequence of laser pulses 
and without individual laser addressing of the ions 
But are these schemes really suitable for quantum com- 
puting with many qubits? Finding reliable ways to scale 
present schemes to many qubits requires simplifications 
of the experimental setup without decreasing the preci- 
sion of gate operations. 

In this paper we present an alternative scheme for 
the realisation of two-qubit gate operations between cold 
trapped ions. The ions are stored inside a linear trap 
(see Figure^ and can be coupled via a common vibra- 
tional mode. Each qubit is, as in Q, obtained from two 
different ground states of the same ion and the system re- 
mains during the whole computation in the ground state 
of the common vibrational mode. Quantum gate oper- 
ations can be implemented within one step by applying 
different laser fields simultaneously. Depending on the 
gate operation, this requires individual addressing of the 
ions by one or two different lasers. Generalising the un- 
derlying idea to the many-ion case might enable us in the 
future to construct efficient schemes for the realisation of 
many-qubit gates and the generation of highly entangled 
many-particle states within one step 

A disadvantage of the proposed scheme is that it allows 
for high gate success rates only in the case of very long 
gate operation times. If the gate is performed fast, so 
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FIG. 1: Schematic view of the experimental setup. Each qubit 
is obtained from two different ground states of the same ion. 
The ions are stored in a linear trap, cooled in the ground state 
of a common vibrational mode and manipulated with the help 
of laser fields individually addressing the ions. 



that its speed becomes comparable to other schemes @1, 
the gate success rate drops for certain initial states below 
80%. However, this disadvantage might be compensated 
by the fact that the scheme is widely robust against pa- 
rameter fluctuations, especially against variations of cer- 
tain Rabi frequencies and detunings. In addition, the 
lasers for the cooling of the common vibrational mode 
are in the following applied continuously which should 
help to decrease decoherence effects due to heating ^ . 
In the setup considered here, there is no need for sympa- 
thetic cooling of the ions with the help of a second ion 
species, as discussed in [^Q- 

The role of cooling in this paper is twofold. On one 
hand, it decreases the sensitivity of the proposed scheme 
with respect to heating. On the other hand, the presence 
of the cooling laser introduces an auxiliary decay channel 
into the system whose presence restricts the time evolu- 
tion of the system onto a small subspace such that the 
implementation of quantum gate operations within one 
step becomes possible. Photon emissions due to cool- 
ing are nevertheless negligible since the system remains 
continuously in a decoherence- free state 10, Jl, 12\ . If 
heating populates the vibrational mode or gate failure 
moves the system out of the decoherence-free subspace. 
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then photons are emitted at a high rate |l3l.ll4l|. This can 
be detected and the computation can be restarted. For 
practical purposes, only the fidelity of gate operations 
must be close to one while the success rate can be signif- 
icantly smaller. As long as the experiment is repeated if 
necessary, quantum computers can still be more efficient 
than conventional computers. 

To achieve this, the scheme utilises an environment- 
induced quantum Zeno effect 0, originating from 
the measurements continuously performed by the cool- 
ing lasers on the common vibrational mode. This is 
possible since ground state cooling has the same effect 
as continuous measurements whether the ions are in the 
n = phonon mode or not. The underlying concept of 
the quantum computing scheme proposed here became 
known in the past as quantum computing using dissipa- 
tion |0| and has especially been used to construct simple 
and precise Ts'l but also very robust quantum com- 
puting schemes for atom-cavity systems. 

Many schemes for quantum computing with trapped 
ions have already been proposed. Some of them re- 
quire cooling of the ions into the ground state of a com- 
mon vibrational mode pj . That this is possible has 
been demonstrated in several experiments 0, 113, |^ . 
Other schemes can be implemented with "hot" ions 
[23. 123*. ^2^. . For example, the proposal by S0renson 
and M0lmer |2^] has been used to entangle up to four ions 
[2^ and to observe a violation of Bell's inequality with 
massive particles jl^- However, most of these schemes 
are very sensitive to level shifts and require precise con- 
trol of the ion-phonon interaction. 

The paper is organised as follows. In the next section, 
we review the basic concept of quantum computing using 
dissipation. Concrete realisations of two-qubit gate oper- 
ations with cold trapped ions are discussed in Section III. 
As examples we consider possible implementations of the 
CNOT gate, a two-qubit phase gate and the SWAP oper- 
ation. Section IV analyses the robustness of the proposed 
quantum computing scheme with respect to fluctuations 
of detunings and other system parameters. Finally, we 
conclude our results in Section V. 



II. QUANTUM COMPUTING USING 
DISSIPATION 

An important feature of the proposed quantum com- 
puting scheme with cold trapped ions is the auxiliary dis- 
sipation channel introduced into the system by continu- 
ous laser cooling. Ground state cooling is a process where 
a detuned laser field transfers the ions into an excited and 
highly unstable atomic state whenever the phonon mode 
becomes populated. During spontaneous emission from 
this level there is a high probability for annihilation of a 
phonon so decreasing the energy in the vibrational mode. 
If the system returns into the ground state with n — 
phonons, the ions no longer see the cooling laser due to 
its detuning to the vibrational side band [a, 0, ^3 and 



the time evolution of the system can be described by a 
Schrodinger equation. 

Although cooling is applied continuously, spontaneous 
emission from the ions is negligible. The reason for 
this is that the state of the system remains during the 
whole computation in a decoherence-free subspace (DFS) 
[TolllllIT^ and no other states become populated during 
gate operations. The DFS of a system is a subspace of 
states whose population does not lead to decoherence. In 
the presence of cooling, the system of cold trapped ions 
possesses a non-trivial DFS with respect to spontaneous 
photon emission. It includes all states with no phonons 
in the vibrational mode while the ions are in a qubit state 
or in a highly entangled auxiliary state. 

To calculate the decoherence-free (DF) states of the 
system and the effective time evolution of the ions dur- 
ing gate operations we use in the following the quantum 
jump approach p^ . This method is equivalent to the 
Monte Carlo wave-function approach 29J and the quan- 
tum trajectory method [30j | and predicts that the initial 
atomic state {ip) evolves under the condition of no photon 
emission as 

\^°{t)) =U,ond{t,0)m . (1) 

Here t/cond(^, 0) is the no-photon time evolution oper- 
ator and corresponds to the non-Hermitian conditional 
Hamiltonian -ffcond [3- For convenience, i?cond has been 
defined such that 

Po{t,ij)^\\U,ond{t,0)\ij)f (2) 

is the probability for no emission in {0,t). 

The basic assumption of the quantum jump approach 
is that the effect of the environment on the state of 
a quantum optical system is the same as the effect of 
rapidly repeated measurements on the free radiation field 
whether a photon has been emitted or not jsj. The 
non-Hermiticity of the conditional Hamiltonian and the 
continuous decrease of the amplitude of the state vector 
\il''^{t)) reflect that the observation of no photons reveals 
information about the system. The longer no photon is 
emitted the more unlikely it becomes that there is exci- 
tation that might cause an emission and the amplitude 
of a state with a spontaneous decay rate decreases expo- 
nentially in time js^ |3^ . 

Using the quantum jump a ppr oach, a state lip) is DF 
if Po{t, V) = 1 for aU times t^n^. Hence, the DFS of a 
system is spanned by the eigenvectors of the conditional 
Hamiltonian with real eigenvalues A^. The eigenvectors 
\Xi) of -ffcond are in general non-orthogonal. It is therefore 
useful to introduce the reciprocal basis vectors jA-') with 
(A-'|Ai) — Sij and to write the conditional Hamiltonian 
as 

J^cond = ^ A, I A,) (A' I . (3) 

i 

Suppose that all non-DF states couple strongly to the en- 
vironment and populating them leads typically to a pho- 
ton emission within a time At. This time At is greater 
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than a certain minimal size which can be determined from 
the quantum jump approach. Provided that the eigenval- 
ues Afc corresponding to non-DF states fulfil the condition 

g-iA.At/ft = , (4) 

the no-photon time evolution operator becomes [34| 



C/cond(At,0) = 



i:|Ai)eDFS 



(5) 



This operator projects every state onto the DFS. The 
action of the environment over a time can therefore 
be interpreted as a measurement whether the system is 
DP or not. The probability for no emission in Ai equals 
the probability to be in a DP state. 

The basic idea of quantum computing using dissipation 
is to utilise the no-photon time evolution jS} for the im- 
plementation of gate operations. As we see below, the 
continuous measurements caused by the environment on 
the quantum system lead to a realm of possibilities to 
induce a DP time evolution between the qubits. The 
reason for this is that the conditional Hamiltonian i?cond 
can easily be varied by an arbitrary but weak enough 
interaction without loosing the restriction of the system 
onto the DPS. As long as the typical time scale of the 
additional terms in the conditional Hamiltonian is much 
longer than At defined by condition Q 17], the result- 
ing effective time evolution is given by © to a very good 
approximation and can be used to generate a unitary op- 
eration between the qubits. 

The restriction of the system onto the DPS can in- 
tuitively be understood with the help of the quantum 
Zeno effect Within At, a weak interaction can 

only transfer population proportional to At out off the 
DPS. During the next measurement the system is found 
in a non-DP state with a probability proportional to At^. 
Otherwise it is projected back onto the DPS by the time 
evolution operator (0). The probability to find the sys- 
tem always in a DP state, i.e. T/ At times if T is the gate 
operation time, goes in the limit of weak interactions, i.e. 
for At/T 0, to one. All transitions out of the DPS are 
therefore strongly inhibited. 

However, the time evolution inside the DPS is not af- 
fected. Using first order perturbation theory with respect 
to the weak interaction, ((Sj and the assumption that the 
system is initially in a DP state, one can show that the 
conditional time evolution in At is the same as time evo- 
lution obtained from the Hamiltonian 



HcS — JPdFS Hcond -^DPS 

using the same approximations. Here 

iPnFS - V I A,) (A, I 



■^DFS - J2 

j:|Ai)eDFS 



(6) 



(7) 



is the projector onto the DPS. The effective Hamiltonian 
-ffoff is used in the following to find the appropriate laser 



configuration for the realisation of certain gate opera- 
tions. 

The reason for the simplicity of the resulting quan- 
tum computing schemes is that the DPS of the cold-ion 
system contains in addition to all ground states, highly 
entangled states. Through populating these states en- 
tanglement between qubits can be created during the ef- 
fective time evolution even if this would not be possible 
using the weak interaction alone. In the next section we 
describe as an example the realisation of universal two- 
qubits gates within one step by applying different laser 
fields simultaneously. 

Deviations from the effective time evolution occur if 
the applied interaction is not weak enough and the sys- 
tem moves in At slightly out of the DPS. We denote 
the (non-unitary) correction to the desired effective time 
evolution in the following by —Acoi-i{T, 0) so that 



f/co„d(T, 0) = [/eff (T, 0) - Acorr(T, 0) 



(8) 



Prom 101 we then find that the gate success rate Po{T, -iji) 
equals 

Po(r,V) - 1 - 2Re(7/.|C/eff(r,0)Acorr(T,0)|7/.) (9) 

while the fidelity of a gate operations under the condition 
of no photon emission is given by 



F(T,V) = 



|(^|t/eff(r, 0);7cond(T, 0)1^)1^ 
Po(T,V) 



(10) 



Using Pquation © it can be shown that the fidelity 
equals one. 



(11) 



in first order in the correction Acorr(T, 0). Dissipation 
can indeed be used to implement very precise quantum 
gates with success rates sufficiently close to one. 

Let us now estimate the probability of finding the re- 
sult of a whole computation assuming that each gate can 
be performed with maximum fidelity but only with a fi- 
nite success rate Pq- The probability of implementing 
an algorithm of N gates faultlessly is Pq and decreases 
exponentially with N . On the other hand, if one always 
knows whether an algorithm has failed or not, the com- 
putation can be repeated until a result is obtained. The 
probability for not having a result after M runs equals 



no result 



= 1 - P 



N 



(12) 



Por large N this is approximately exp(— MPj^). Many 
repetitions might be necessary to implement a computa- 
tion. However, for smaller numbers of N and if Pq is suf- 
ficiently close to unity, the failure probability is already 
nearly negligible for M N . Por example, if Pq = 95% 
and an algorithm with A'^ = 50 gates is performed, then 
repeating the computation 50 times yields already a suc- 
cess rate above 98%. 
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FIG. 2: Level scheme of the two ions involved in the gate 
operation. Each qubit is obtained from the ground states 
|0) and |1) of one ion. In addition a metastable state |2), a 
rapidly decaying state |3) with decay rate Fa and two strong 
laser fields with coupling strength gj = ^rjjQj and detuning 
V are required. 

III. SINGLE LASER PULSE GATES 

In this section we propose concrete realisations of two- 
qubit quantum gates for cold trapped ions. Each qubit 
is obtained from two different atomic ground states |0) 
and 1 1) of the same ion. In addition, a metastable state 
1 2) and a rapidly decaying level 3 are required, as shown 
in Figure 121 First, the ions have to be cooled into the 
ground state of a common vibrational mode. Two strong 
laser fields detuned by the frequency of a common vi- 
brational mode should be applied. The laser field cou- 
pling to the 1-2 transition establishes an "interaction" 
between the two qubits involved in the gate operation. 
The laser field driving the 1-3 transition represents the 
laser cooling setup and can be replaced by any other laser 
cooling configuration without changing the effective time 
evolution of the system. 

In the following, we denote the spontaneous decay rate 
of level 3 by while b and are the annihilation and 
creation operator of a phonon in the common vibrational 
mode. The coupling constant of this mode to the atomic 
1-j transition equals gj = ^-qj^lj where ftj is the Rabi 
frequency of the applied laser field and rjj is the Lamb 
Dicke parameter depending on the characteristics of the 
ion trap. Proceeding as in one finds that the con- 
ditional Hamiltonian within the dipole and the rotating 
wave approximation and in the interaction picture with 
respect to the free Hamiltonian equals 

2 

Hcond = ^i/l[52|l),(2|&U53|l)»(3|6^-h.C.] 

i=l 

2 

-^iftr3|3).(3| . (13) 

i=l 

Here the Lamb-Dicke regime and the condition i' ^ il,2, 
^3 has been assumed, as in [J. 

The DF states of the two ions are spanned by the eigen- 
states of the conditional Hamiltonian iJcond with no pop- 
ulation in the unstable state |3). The DFS of the system 
therefore contains only superpositions of states with the 



ions either in the state |00) combined with an arbitrary 
state of the vibrational mode or both ions in |01), |10), 
1 11) or in the antisymmetric state 

|a)^^[|12)-|21)] (14) 

while the vibrational mode is not populated. The DF 
states correspond indeed to the eigenvectors of iJcond 
with real eigenvalues. Here these are all zero eigenvalues 
and the system does not evolve as long as no additional 
interaction is applied. 

To realise gate operations between the two qubits weak 
laser fields are required in addition to the strong lasers 
shown in Figure |21 Let us denote the Rabi frequency of 
the laser with respect to the j-2 transition in ion i by 
ri^*-* and assume that 

nf « <?2, 53 and Fa . (15) 
The laser Hamiltonian equals then 

= E E 5^-'^b'>^(2|+h.c. (16) 

i=1^2 j=0,l 

If the system is restricted onto the DFS, as predicted in 
the previous section, then the time evolution of the ions 
can be predicted with the effective Hamiltonian (jSJ which 
yields 

-n['^)\u)]{a\ + h.c. (17) 

The operation time T should be chosen such that all pop- 
ulation returns again into the qubit states at the end of 
each gate operation. 

To analyse the no-photon time evolution of the ions 
in more detail and to determine the optimal parameter 
regime for quantum gate implementations the notation 

oo 3 

= E E ^'Mji^^^i) (18) 

is introduced. Here Cn,ij is the amplitude of the un- 
normalised state with respect to the state with n 
phonons in the vibrational mode while the ions are in 
\ij). Under the condition that 

fl!f « 32, 53 and Fa , (19) 

the coefficients Cn,ij of all non-DF states change on a 
much faster time scale than the coefficients of the DF 
states. They can therefore be eliminated adiabatically 
by setting their derivatives equal to zero. Using lfn|l . 
(|16|l and the Schrodinger equation for the no-photon time 
evolution we obtain 
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9l 



45|r3 



1 



9l 



gl 



1 



(21) 



As long as the parameters h\ and fc2 muhiphed with the 
Rabi frequencies oj'-* in the effective HamiUonian p7|l 
are much smaller than one, the system behaves indeed as 
predicted by the quantum Zeno effect. 

In the following subsections we discuss in detail three 
concrete quantum gate implementations, starting with 
the CNOT gate. This gate, as well as the phase gate, 
constitutes together with the single-qubit rotation a uni- 
versal set of quantum gates [33 ■ Another gate, which 
as we will see can easily be performed with a very high 
fidelity and success rate is the SWAP operation. 



A. CNOT gate 

If ion 1 contains the target qubit and ion 2 provides 
the control qubit, then the CNOT gate is given by the 
time evolution operator 

f/gate = |00)(00| + |01)(01| + |10)(11| + |11)(10| . (22) 

The easiest way to realise this gate is when one laser cou- 
ples with the (real) Rabi frequency 17 to the 1-2 transition 
of ion 1 and another one with the same Rabi frequency 
to the 0-2 transition of ion 2, i.e. 



d-^^ = d-^^ = n and n};-' = ni-"' = o 



(1) _ o(2) _ 



(23) 



This parameter choice corresponds to the effective Hamil- 
tonian 

Hes = ^;ir![|10>-|ll>] (a|+h.c. (24) 

If the duration of the laser pulse equals 

T = 2n/n, (25) 

the resulting time evolution is exactly the desired CNOT 
operation. 

From (|20|l and first order perturbation theory we find 
that the correction of the no-photon time evolution to 



the effective time evolution of the DF states equals in 
first order in kiVl 

Acorr(T,0) = 2r2^fciT|00)(00| 

+ \n^{7ki + k2)T [|10>(11| + h.c] 
+ ir!2(5fci - k2)T [|10)(10| + |11)(11|] 
-2V2iflki [\10){a\ + |ll)(a| +h.c.] 
+ [l-^n''{ki+k2)T]\a){a\ . (26) 

This leads to the no-photon probability and gate success 
rate 

Po(r,V) = l-4r!2fciT|co,oo|' 

-n2(5fc^_fc2)rRe[co,iocS,n] 
-^n''i7ki + k2)T [|co,ioP + |co,iiP] 
-4V2nki Im[co,ac5,io + co.qCq (27) 

which is very close to one as long as fciT and k2T are 
small. Using H26|l one can calculate the gate fidelity under 
the condition of no photons. As expected, it equals one 
in first order in kiQ. 

Let us now determine the parameters that maximise 
the minimal gate success rate PoiT^ip) on average, this 
means independent from the initial qubit state. From (j^J 
we find that 



Po{T, ip) > min 



(28) 



if fii are the eigenvalues of the conditional time evolution 
operator p6() . These are always smaller but close to one. 
Calculating them we find that up to first order in ki one 
has 



Po{T,i;) > 1 -n^ina^{{ki + /cs)^, 6kiT} 



(29) 



Optimising this we find that for a fixed ratio of Q / g2 the 
CNOT gate works most reliable if 



53opt = V2g2 and Fsopt = 2\/37 g2 



(30) 



For these parameters one has, for all qubit states 

^Oopt(T, V) > 1 - 37rn/{VS7g2) . (31) 

The success rate of the gate operation assumes its mini- 
mum if is a superposition of the states |10) and |11). 



6 



0.9 

Po(T,V) 

0.8 







(a) - 






(b) 






(c) 






(d) 













0.05 0.1 0.15 0.2 

Q/g2 



F(T,v) 




0.988 - 



0.984 I 1 1 1 J 

0.05 0.1 0.15 0.2 

a/g2 

FIG. 3: Success rate and fidelity under the condition of no 
piioton emission of a single CNOT gate as a function of the 
Rabi frequency il for Fa = Fsopt , ~ gsopt and for the initial 
qubit states |00) (a), |10) (b), |11) (c) and [|10) - |11)]/^ 
(d). The gate success rate and fidelity are always maximal if 
the ions are initially prepared in |01). 
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FIG. 4: Fidelity for a single CNOT gate under the condition 
of no photon emission as a function of the spontaneous decay 
rate Vz/g2 for = Fs, Q = 0.01 32 and for the initial qubit 
states |00) (a), |10) (b) and [|10) - |11)]/^ (c). If the ions 
are initially prepared in |01), the gate success rate equals one. 
For — 1 10) the fidelity is about the same as in graph (c). 

cial for the scheme to work! For small damping rates the 
mechanism which restricts the time evolution of the sys- 
tem onto the computational subspace fails and the min- 
immii gate fidelity is well below 50%. If performed with 
ions with a larger spontaneous decay rate T3, the scheme 
works precisely - although this parameter should also not 
be too large since this would introduce an additional fast 
time scale into the system and the assumptions made to 
derive the differential equations H2(J|) would no longer be 
valid. 



For all other initial qubit states the success rate is even 
closer or equal to one. 

Figure |2| results from a numerical solution of the no- 
photon time evolution given by the sum of the Hamil- 
tonians (|13|l and 116|l . The success rate Pq is in good 
agreement with the theoretical predictions H31|) . For very 
small Rabi frequencies, the gate success rate and fidelity 
is for all initial states close to one. For larger values of Jl, 
the gate success rate decreases and can for V, = 0.2 52 be 
as low as 73%. The gate fidelity is in this case still above 
98.4%. The smallest gate success rate is found when the 
atoms are initially in |10), or in a superposition of 
these two states. Success rates Pq > 90% are achieved as 
long as < 0.07 32- The fidelity F{T,^) is in this case 
larger than 99.8% and the gate duration time T is about 
90/32. 

Finally, we would like to comment on the role of the 
spontaneous decay rate of level 3 in the scheme. Figure^ 
shows the fidelity of a single CNOT gate as a function of 
T3 and for — T3. In the chosen parameter regime, the 
effective damping rate of unwanted population in non- 
DF states can be shown to equal g^/T^ = Fa to a very 
good approximation and the effective decay rate of non- 
DF states increases linearly in T3. Figure^confirms that 
the presence of the auxiliary dissipation channels is cru- 



B. Phase gate 

Another two-qubit quantum gate, whose realisation is 
even simpler than the realisation of the CNOT gate, is 
the phase gate that changes the sign of the amplitude of 
the qubit state |01) but leaves all other basis states un- 
changed. It requires only one weak laser pulse in addition 
to the laser excitation shown in Figure |21 3 ■ This pulse 
addresses the 0-2 transition in atom 1 and one should 
choose 

n^^^ = n and nl^^ = = of ^ = o . (32) 

From (|17|l we see that this choice of parameters corre- 
sponds to the effective Hamiltonian 

H,s = ^hn\01){a\+h.c. (33) 
and the weak laser field should be turned off after a time 

T = 2\/27r/l] . (34) 

To change the sign of the qubit state |10) one can proceed 
in exactly the same way but with fig^^ = f2 and f^Q^-* = 
f]«=f]f)=0. 
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FIG. 5: Success rate of a single phase gate and fidelity under 
the condition of no photon emission as a function of the Rabi 
frequency fl for the same parameters Fa and gs as in Figure 
|3| The initial states of the ions are |00) (a), |01) (b) and 
[|00) + |01)]/V2 (c). 



Using the differential equations H20I) and first or- 
der perturbation theory we find that the correction 
^corr(T', 0) to the effective time evolution equals in first 
order kifl 

Acorr(T,0) = %/27rfcirj [4|00)(00| - |01)(01|] .(35) 



C. SWAP operation 

Another quantum gate that can be implemented rel- 
atively easy is the SWAP operation. It transforms the 
state |01) into |10) and vice versa but leaves |00) and 
1 11) unchanged. The realisation of this gate does not 
require individual laser addressing of the ions and one 
should choose 

n^^^ ^ nl^^ = n and = f^f ^ = o . (3?) 

For these Rabi frequencies the effective Hamiltonian H17() 
becomes 

H,ff = ^hn[-\01) + \10)]{a\+h.c. (38) 

and the duration of the weak laser pulse should equal 

T = 2TT/n . (39) 

The SWAP gate can be very useful. It exchanges the 
states of two qubits without that the corresponding ions 
have to swap their places physically. 

Proceeding as in the previous subsection, the correc- 
tion to the desired effective time evolution reads 

Acorr(r,0) = 27rfcil7 [4|00)(00| + |01)(01| + |10)(10| 
+ |01>(10| + |10)(01|] ® |0)vib(0| (40) 

and the gate success rate becomes 

Po(T,V) - l-4^fciri[4|co^oo|' + |co.oi|' + |co,io|' 
+2Re(co,oi4io)] ■ (41) 

Again, the fidelity under the condition of no photon emis- 
sion can be shown to equal one in first order fcifi. A 
comparison with H27|) shows that the success rate of each 
SWAP gate is of similar size as for the phase gate dis- 
cussed in the previous subsection. 



This leads to the no-photon probability 

Fo(r,V) = l-2V27rfcil][4|co,oo|' + |co,oin (36) 

which is very close to unity as long as fcifi ^ 1. Calculat- 
ing the fidelity under the condition of no photon emission 
(|1U|1 one finds again that it does not differ from one in 
first order. 

Figure results from a numerical solution of the no- 
photon time evolution of the system for the same param- 
eters as in FigureUI The success rate and the conditional 
fidelity of a single phase gate are of similar size as in the 
previous subsection and the scheme works very well for 
small Rabi frequencies fJ. As initial states we considered 
only the the qubit states |00) and |01) since they are the 
only ones affected by the laser field. If the ions are ini- 
tially in 1 10) or |11), the fidelity and success rates of the 
phase gate are exactly one. 



D. Dissipation-assisted quantum gates 

A closer look at the differential equations 1)20(1 reveals 
that for fixed gate operation times relatively small cor- 
rections to the effective time evolution (fT7j) are obtained 
if 

f]W _ ^7(2) ^ Q ^ ^ ^^2) 

Then all terms proportional to ^2 vanish in 1(20(1 and the 
terms proportional to fci become relatively small. The 
corresponding gates are relatively easy to realise because 
they require only one weak laser field applied to the 0-2 
transition in each atom. The effective Hamiltonian is in 
this case given by 

H,s = ^n[-ni^^\oi) + nl,^^\m)]{a\ + h.c. 

(43) 
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n/g2 

FIG. 6: Fidelity for a single phase gate under the condition 
of no photon emission as a function of the Rabi frequency 
Q if the ions are initially prepared in the state |00) and for 
(/3 = Fs = (a), for 33 = (?2 and F3 — W g2 (b) and for 
ff3 = 92 and r3 = (c). The curves give a lower bound for 
the fidelity of the prepared state independent from the initial 
state. 



and can be realised with fl\ — fl\ — 0. The above 
described phase gate and SWAP operation are examples 
for quantum gates that can be implemented in this way 
and they are the only gates for whom the presence of the 
auxiliary dissipation channel is not crucial for the scheme 
to work. 

Remarkable about the effective time evolution (|43|l is 
that cooling plays indeed a different role for its realisa- 
tion than in the case of the CNOT gate which only works 
because of dissipation. Suppose the quantum Zeno effect 
provides the basic mechanisms that restricts the system 
onto the DFS. Then one would expect that the quantum 
gate works the better if one increases the value of /T^ 
since this number relates directly to the spontaneous de- 
cay rate of non-DF states. Instead it can be shown that 
the phase gate and the SWAP operation work the better 
the smaller /Fa since this implies smaller values for fci . 
Both gates can even be implemented in the absence of 
cooling for g^ — T3 — 4]. 

The mechanism which yields the effective Hamiltonian 
H43|l in this subsection is adiabaticity resulting from the 
presence of two different times scales in the system. One 

(i) 

is given by the weak Rabi frequencies f2g and the other 
one by the coupling constant g2. To guarantee a suc- 
cess rate close to one, quantum gates based on adiabatic 
passages are as slow as gates based on a quantum Zeno 
effect. Advantages arise when the system is operated 
outside the adiabatic regime with gate operation times 
of about the same size as 1/(72- In this regime, cool- 
ing has the same effect as error detection measurements. 
Under the condition of no photon emission, the system 
behaves as predicted by adiabaticity. The presence of the 
auxiliary dissipation channel exponentially damps away 
the population that accumulates due to non-adiabaticity 
in unwanted states. The resulting process can be called a 
dissipation-assisted adiabatic passage The price one 
pays for the speed up of gate operations is a decrease of 



the gate success rate Pq. 

Figure shows the fidelity of a single phase gate under 
the condition of no photon emission for a wider range of 
Rabi frequencies fl than assumed in FigurcsElandEl The 
fidelity of the finally obtained state is for 53 = F3 = 
very close to one as long as fl is much smaller than 52 
(see FigureEIa)). Note that the same fidelity close to one 
can be achieved for higher Rabi frequencies, and therefore 
also for shorter gate operation times, in the presence of 
dissipation (see FigureEfb)). Choosing g^ — g2 and F3 — 

improves the precision of the performed gate operation 
even further, corresponds to the gate success rate Pq = 

1 but does not provide cooling of the ions during gate 
operations (see Figure Efc)). 



IV. ROBUSTNESS AGAINST PARAMETER 
FLUCTUATIONS 

In this Section we discuss in detail the robustness of 
the proposed quantum computing scheme with respect 
to parameter fluctuations and analyse the effect of varia- 
tions of Rabi frequencies and the presence of detunings. 
The independence of the scheme from the precise value 
of most experimental parameters is a big advantage of 
the proposed quantum computing scheme and applies to 
quantum computing schemes using dissipation in general. 



A. Effect of detunings 

In the following we denote the detuning of the laser 
field that establishes the coupling of the ions to the com- 
mon vibrational mode by A2. The reason for the inde- 
pendence of the no-photon time evolution of the system 
from this parameter arises from the different role of the 
ion-phonon interaction in the present scheme compared 
to other quantum computing schemes with cold trapped 
ions. Here the interaction does not cause a transfer of 
excitation into a common vibrational mode. Instead, 
its presence introduces a fast time scale into the system 
which prevents transitions into unwanted states. During 
gate operations, the phonon mode becomes only virtually 
populated. 

To visualise the role of the vibrational mode in the 
scheme it is useful to have a look at all levels and transi- 
tions involved in the time evolution of the system. As an 
example we survey in the following the phase gate which 
has been introduced in Section IIII Bl Its level scheme 
is shown in Figure |7| and is much simpler than the level 
configurations for the CNOT gate and the SWAP opera- 
tion. Still, the independence of detunings described here 
applies to single laser pulse gates with cold trapped ions 
in general and can be shown for all other gates in the 
same way. 

If initially only qubit states are populated, then the 
time evolution of the system during a phase gate imple- 
mentation involves only the eleven levels shown in Figure 
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FIG. 7: Level scheme showing all transitions and level in- 
volved in the realisation of the two-qubit phase gate. 



[7| Let us define the interaction-free Hamiltonian Hq as 

Ho = -ffions + Hp\ 

lonons 

?iA2[|l,10)(l,10| 
+ |1, 11) (1, 111 + |0, 30) (0, 30| + |0, 513) (0, sial ] , 

(44) 

where A2 is the detuning of the strong laser field driving 
the 1-2 transition of each ion and 



\s,,)^^[\lj) + \jl)] 



(45) 



In the interaction picture with respect to this Hamilto- 
nian, the Hamiltonian describing the no-photon time evo- 
lution of the system becomes 

Hcond = i/i52 [|0,20)(l,10| + x/2|0,si2)(l,ll|] 
+i%3 [ |1, 10)(0, 30| + \/2 |1, 11)(0, S13I ] 

[ |o, 00) (0, 20K ^ |o, 01) (0, S12I 

+ ^|0,01)(0,a|] +h.c. + nA2 [|1,10)(1,10| 

-Mi,ii)(i,ii|] -?i(ir3- A2) [|o,30)(o,30| 

+ |0,si3)(0,si3|] . (46) 

One can easily check that the effective Hamiltonian © 
of the system does not depend on the detuning A2. In 
a first approximation, the time evolution of the system 
does not change, even if the detuning A2 is about the 
same size as 32- 

This is in good agreement with the numerical results 
shown in Figure|Slwhich has been obtained by solving the 
no- photon time evolution of the system with the Hamilto- 
nian (|44|l . As expected, the presence of the detuning A2 
has only a small effect on the success rate of the scheme. 
The fidelity of the gate operation under the condition of 
no photon emission decreases less than one percent even 
if A2 becomes as big as the ion-phonon coupling constant 
g2. However, the amplitude of the finally prepared state 
collects a phase factor in the presence of detuning. De- 
pending on the initial qubit state, this phase can differ 
up to 12° from the desired phase (see Figure El . 

Analogously, it can be shown that detunings in the 
cooling laser setup affect the effective time evolution of 
the system only slightly. Detunings of the weak laser 
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FIG. 8: Success rate of a single phase gate and fidelity under 
the condition of no photon emission as a function of the de- 
tuning A2 of the ion-coupling laser for the same parameters 
and initial states as in Figure |^ The last graph gives the size 
of the phase error Aip of the amplitude of the final state in 
degrees. 



field with Rabi frequency Q have to be taken into account 
because they lead to an additional term in the effective 
Hamiltonian 



B. Fluctuations of Rabi frequencies 

Let us now assume again zero detunings and discuss 
the effect of fluctuations of Rabi frequencies on the de- 
sired time evolution of the system. The analysis in Sec- 
tion ^O] showed already that the performance of the pro- 
posed quantum gate operations does not depend on the 
size of the ion-phonon coupling constants (72 and (73 and 
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the decay rate T3 as long as these rates are big enough. In 
the formahsm this is reflected by the fact that the effec- 
tive Hamiltonians (|24|l . and (|38|l depend only on fl. 
In this subsection we show that the scheme is also widely 
robust against fluctuations of the weak Rabi frequency 

The reason for this robustness is that the effective 
Hamiltonian can be written as HcS = where M is a 
time independent operator. This is the case for all three 
quantum gates discussed in the previous Section. Let us 
assume that the Rabi frequency fl = il{t). The effective 
time evolution operator of the system is then given by 



C/eff(T,0) = exp 



dtn{t)M 



(47) 



and depends only on the integral dtn{t). High fi- 
delities require therefore only precise control of this time 
integral. With respect to fluctuations of system parame- 
ters, quantum computing schemes using dissipation share 
some features with schemes that employ holonomies or 
geometrical time evolutions j^Bl or similar dynamical pro- 
cesses to implement quantum gate operations. 

Finally, we would like to comment on the effect of 
shaped laser pulses on the quantum gate performance. 
As an example we consider again the phase gate and as- 
sume 



n{t) = 2^0 sin^ 



while the gate operation time should be chosen as 



(48) 



(49) 



The integral dtfl{t) is then the same as in Subsection 
IIII Bl and the Rabi frequency equals zero at the end 
of the gate operation. Motivation for assuming a time 
dependence as in (|48(l is to minimise the corrections to 
the desired time evolution and to increase the fidelity 
of the finally prepared state without increasing the gate 
operation time. 

Figure |51 results from a numerical simulation of the 
no-photon time evolution of the system and shows that 
F{T, ip) can be improved even further. The gate fidelity is 
now for many initial states closer to one and for all initial 
states above 99.2%, even if flo = 0.2 g2. A comparison of 
Figure El with Figure El reveals that the price one has to 
pay for this improvement of the precision of the scheme 
is a reduction of the gate success rate Pa{T, for certain 
initial qubit states by a few percent. (Note that the gate 
operation times T are in both Figures the same when f2 
and r^o a-rc the same.) Similar results can be derived for 
other single-laser pulse quantum gates with cold trapped 
ions. 



V. CONCLUSIONS 

We discussed the possibility to implement precise 
quantum gates between cold trapped ions in the presence 



Po(T,¥) 0.8 



F(T,xt/) 




0.992 



FIG. 9: Success rate of a single phase gate and fidelity under 
the condition of no photon emission as a function of the Rabi 
frequency 57o (defined in 1481 1 ) for the same parameters and 
initial states as in Figure |S] 



of cooling of a common vibrational mode within one step. 
That the proposed quantum computing scheme is feasible 
with present technology has recently been demonstrated 
in an experiment in Innsbruck 0. The setup for the 
realisation of the Cirac-ZoUer CNOT gate with Calcium 
ions reported in is very similar to the setup needed for 
the implementation of the two-qubit phase gate and the 
SWAP operation presented in this paper. Both quantum 
computing schemes require cooling of the ions to the vi- 
brational ground state and a strong laser field to establish 
a coupling between the ions. In addition only one laser 
pulse is needed. To implement the phase gate, this laser 
has to excite one of the ions resonantly. 

Because of its simplicity and robustness, the proposed 
scheme might help to increase the number of qubits in 
present quantum computing experiments. One of its ad- 
vantages is that it offers a great variety of ways to imple- 
ment gate operations, which can all be performed within 
one step. Laser cooling is not only possible between com- 
putations but should be applied continuously. An addi- 
tional ion species for sympathetic cooling, as proposed in 
0, is not required. Although the system remains con- 
tinuously in the ground state of a common vibrational 
mode, the presence of the cooling laser setup decreases 
the probability for heating even if no photons are actually 
scattered from the ions [6| . 

Another major benefit of the proposed quantum com- 
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puting scheme is its independence from fluctuations of 
most experimental parameters. To guarantee a high pre- 
cision of gate operations only one of the applied laser 
fields has to be controlled relatively well. Its Rabi fre- 
quency integrated over the gate operation time has to 
equal a certain value. All other Rabi frequencies, Q2 and 
r^s, do not enter the effective time evolution of the sys- 
tem and the scheme does not depend on their concrete 
value as long as they are big enough. The proposed setup 
can also tolerate a detuning of the laser field which estab- 
lishes the coupling of the ions to the vibrational mode, 
even if this detuning is of the same size as the ion-phonon 
coupling constant (?2- 

Error detection is already included in the scheme Q. 
If a gate operation fails or heating leads to the popula- 
tion of higher phonon modes, then photons are emitted 
at a high rate. This can easily be detected and the whole 
computation can be restarted if necessary. In this way, 
the system is protected against certain kinds of errors. 
As long as the gate success rate is sufficiently close to 
one, quantum computing can still be more efficient than 



performing the same computation on a conventional com- 
puter. 

A disadvantage of the presented quantum computing 
scheme is that its gate operation times are longer than 
in many other setups. For example, the implemen- 
tation of a single two-qubit phase gate requires the 
time T = 2\/2tt/Q where fl can be as large as 0.1 
times the ion-phonon coupling constant 52- Under 
the condition of no photon emission, this guarantees a 
fidelity above 99.6 % and a gate success rate above 85 %, 
independent of the initial state of the ions (see Figure 
Isj. Higher fidelities and success rates are obtained for 
smaller Rabi frequencies and longer gate operation times. 
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